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Marine invertebrates interact with various microorganisms ranging from pathogens to symbionts. One-to-one 
symbiosis between a single microbial species and a single host animal has served as a model for the study of host- 
microbe interactions. In addition, increasing attention has recently been focused on the complex symbiotic associations, 
e.g., associations between sponges and their symbionts, due to their biotechnological potential; however, relatively 
little is known about the microbial diversity associated with members of the phylum Echinodermata. Here, for the 
first time, we investigated microbial communities associated with a commercially important holothurian species, 
Apostichopus japonicus, using culture-dependent and -independent methods. Diverse and abundant heterotrophs, mostly 
Gammaproteobacteria members, were cultured semi-quantitatively. Using the cloning and sequencing technique, 
different microbial communities were found in different holothurian tissues. In the holothurian coelomic fluid, potentially 
metabolically active and phylogenetically unique members of Epsilonproteobacteria and Rickettsiales were discovered. 
This study suggests that coelomic fluids of marine invertebrates, at least those inhabiting intertidal areas where physical 
and chemical conditions fluctuate, provide microbes with unique and stable habitats. 
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MICROBES AND 
ENVIRONMENTS 



Marine invertebrates harbor a wide array of ecto- and 
endo-symbiotic microbes. One-to-one symbiosis between a 
single microbial species and a single host animal has served 
as a model for the study of host-microbe interactions (7, 22, 
37, 38, 39). Recently, increasing attention has been focused 
on complex symbiotic associations between a single host 
species and multiple microbial species/populations, e.g., 
sponge symbiosis due to their biotechnological potential (57, 
58). Many marine invertebrates, particularly sessile to nearly 
sessile, are rich sources of valuable metabolites (9, 14). An 
individual sponge harbors dozens of different microbial 
species, in which bacteria comprise up to 40% of total sponge 
volume (59). Sponge-dwelling bacteria or the sponge itself 
produce a variety of bioactive compounds showing antiviral, 
antitumor, and antimalarial activities (45, 58). From most 
phyla of marine invertebrates, host-associated microbial 
communities have been characterized; however, relatively 
little is known about microbial diversity within and on 
Echinodermata members such as sea stars, sea urchins, and 
sea cucumbers (19). 

Holothurians are a group of marine invertebrates belonging 
to Echinodermata, and are distributed widely throughout 
the world from intertidal zones to the deep sea. Some 
holothurian species have a high market value and thus are 
heavily overexploited (18). Holothurians are model organ- 
isms for studying processes related to organ regeneration and 
organogenesis since they regenerate their digestive tract 
within a few weeks after self-evisceration (12, 42). Their 
body wall and coelomic fluid contain bioactive compounds 
with antibacterial, antifungal, anticoagulant, and antitumor 
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activities (8, 17, 29). Although some microbes, e.g., members 
of the genera Pseudoalteromonas and Salegentibacter, have 
been isolated in pure cultures from holothurians (1, 23, 34, 
48), little is known about the associated microbial diversity. 

In this study, we investigated microbial communities 
associated with a commercially important holothurian, 
Apostichopus japonicus {Echinodermata: Holothuroidea). 
We quantified culturable microorganisms in major micro- 
bial habitats in the holothurian body, including the body 
surface, digestive tract, and coelomic fluid. In addition, 
microbial diversity and niche segregation were comprehen- 
sively and quantitatively evaluated using culture-independent 
approaches. 

Materials and Methods 

Sampling and processing 

Individuals of natural holothurians, A. japonicus (Echinodermata: 
Holothuroidea), were collected from the Funka Bay and Ainuma 
fishing port (42°08'N, 140°06'E), Hokkaido, Japan. Individuals 
(15.0-22.1 cm in length, 195-315.5 g) were rinsed with sterile 75% 
ASW (NaCl, 22.5 g; KC1, 0.5 g; MgS0 4 -7H 2 0, 4.0 g; MgCl 2 -6H 2 0, 
8.1 g; CaS04-2H 2 0, 1.0 g per liter) and were aseptically dissected 
into the small intestine, large intestine, body surface, and coelomic 
fluid. Briefly, outer surfaces of the body wall (12 cm 2 per individual) 
were sampled using sterile cotton-tipped swabs (Men-tip; J.C.B, 
Industry, Tokyo, Japan). After wiping the body surface with 70% 
ethanol, the body wall was cut with a disposable scalpel blade to 
collect the coelomic fluid using a sterile syringe. Finally, small and 
large intestines were sampled. The cotton-tipped swabs and digestive 
tracts were homogenized in sterile 75% artificial seawater (ASW) 
using a stomacher (Pro-media, SH-IIM; ELMEX, Tokyo, Japan) 
(49). Samples were serially diluted in sterile 75% ASW for the 
cultivation test. 
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Table 1. Primers used in this study 



Primer 


Target 


Sequence (5'-3') 


Position (E. coli) 


Reference 


24F 


Bacteria 


AGAGTTTGATCCTGGCT 


8 to 24 


(50) 


1509R 


Bacteria 


GGCTACCTTGTTACGACT 


1,510 to 1,493 


(50) 


1540R 


Bacteria 


AAGGAGGTGATCCAGCCGCA 


1,541 to 1,522 


(50) 


Arch21F 


Archaea 


TTCCGGTTGATCCYGCCGGA 


7 to 26 


(3) 


Arch915R 


Archaea 


GTGCTCCCCCGCCAATTCCT 


935 to 915 


(3) 


920R 


Bacteria 


CCCCGTCAATTCCTTTGAGT 


928 to 909 


(50) 



Cell counting and cultivation test 

Total cell counts were obtained by direct cell counting of 
the glutaraldehyde-fixed cells with DAPI (4',6-diamidino-2- 
phenylindole) using epifluorescence microscopy (47). The abun- 
dance and diversity of culturable heterotrophs were estimated as 
previously described (56). After 7 days' cultivation under 20°C, 
about thirty colonies were randomly picked up from plates. Isolates 
were physiologically characterized as previously described (49). In 
addition, 16S rRNA gene sequences of isolates were determined 
using the primers 24F and 1509R (Table 1). The PCR conditions 
were as follows: 30 cycles of 94°C for 1 min, 55°C for 1 min, and 
72°C for 1.5 min. GoTaq Green Master mix (Promega, Madison, 
WI, USA) was used for the PCR reactions. 

16S rRNA gene clone library 

Microbial diversity was assessed for holothurians collected from 
the Ainuma fishing port. Microbial DNA was directly extracted 
from the microbial community in each homogenate sample or 
coelomic fluid with a Wizard Genomic DNA Purification Kit 
(Promega). DNA extracted from seawater surrounding the holothu- 
rian community was used as a positive control. To reduce PCR bias 
for the assessment of bacterial diversity, two primer sets, 24F-1 509R 
and 24F-1540R, were used (Table 1). The PCR conditions were as 
follows: 35 cycles of 94°C for 1 min, 50°C for 1 min, and 72°C 
for 1 .5 min. Each amplicon was excised and purified (Wizard SV 
Gel and PCR Clean-up System; Promega), and then ligated into 
the TOPO TA cloning vector (Invitrogen, Carlsbad, CA, USA) 
according to the instructions. Ligation products were transformed 
into Escherichia coli One Shot TOP 10 cells (Invitrogen). Clones 
were amplified by PCR with vector-specific primers. The 920R 
primer was used for partial sequencing of the insert to determine 
the phylogenetic clone type (phylotype). Clones with >94% 
similarity were assigned to the same phylotype. Approximately 800 
bp of each representative rRNA gene clone sequence was determined 
for both strands. To estimate the representation of the phylotypes, 
coverage was calculated by Good's equation (13) with the formula, 
(l-[«;/A'])xl00, where «; is the number of single-occurrence 
phylotypes within a library and N is the number of clones examined. 
The bacterial community structures were compared by cluster 
analysis based on the clonal frequency of each representative 
phylotype. The square distance was determined by the Ward method 
(41). 



Construction of phylogenetic tree 

The 16S rRNA gene sequences of representative clones and 
isolates were aligned with ARB software (25). Alignments were 
manually verified with known secondary structure regions. Phylo- 
genetic analyses were restricted to nucleotide positions that could 
be unambiguously aligned. Phylogenetic trees were generated by a 
distance method using PAUP* 4.0b (54) and ARB. Distances were 
estimated with the Jukes-Cantor correction. Bootstrap analyses with 
100 trial replications were used to obtain confidence estimates for 
the tree topologies. 

FISH analysis 

Fluorescent in-situ hybridization was performed as described 
elsewhere (3, 21, 40, 51). In brief, cells were hybridized with the 
Cy3-labeled probes (Table 2) for 4.5 h at 46°C. The percentage of 
fluorescently-labeled cells to DAPI-stained cells was determined 
using epifluorescence microscopy. 

Nucleotide sequence accession numbers 

The 16S rRNA gene sequences of the representative isolates and 
clones obtained in this study have been deposited in DDBJ/EMBL/ 
GenBank under Accession No. AB550432 to AB550557. 

Results and Discussion 

Total cell counts and cultivation test 

Total cell densities were estimated for the holothurian body 
surface, large and small intestines, and coelomic fluid from 
two individuals (Table 3). Unexpectedly, many prokaryotic 
cells were observed even in the coelomic fluid, which varied 
in the range of 4.8><10 5 to 9.9><10 5 cells mU 1 . Cells observed 
in the coelomic fluid were mostly twisted rods (Fig. SI). The 
culturable heterotrophic populations generally accounted for 
less than 5% of the total cell counts (Table 3). In the sample 
Li-2 (the large intestine from the Ainuma individual), the 
population of culturable heterotrophs was probably overes- 
timated, potentially due to the uneven distribution of cells in 
the large intestine homogenate. 



Table 2. Oligonucleotide probes and competitors used in this study 



Probe or competitor 


Specificity 


Sequence (5'-3') 


Target 
RNA 


Position (E. coli) 


FA 

(%) 


Reference 


Probe ARCH9 15 


Archaea 


GTGCTCCCCCGCCAATTCCT 


16S 


915-935 


20 


(53) 


Probe EUB338 


Bacteria 


GCTGCCTCCCGTAGGAGT 


16S 


338-355 


20 


(2) 


Probe ALF968 


Alphaproteobacteria 


GGTAAGGTTCTGCGCGTT 


16S 


968-986 


20 


(35) 


Probe GAM42a 


Gammaproteobacteria 


GCCTTCCCACATCGTTT 


23S 


1,027-1,043 


35 


(26) 


Probe BET42a 


Betaproteobacteria 


GCCTTCCCACTTCGTTT 


23S 


1,027-1,043 


35 


(26) 


Probe EP402-423 


Epsilonproteobacteria 


GAAAKGYGTCATCCTCCACG 


16S 


402-423 


30 


(55) 


Probe CFB560 


Bacteroidetes 


WCCCTTTAAACCCART 


16S 


560-575 


30 


(43) 


Probe PLA886 


Planctomycetes 


GCCTTGCGACCATACTCCC 


16S 


886-904 


35 


(36) 


Competitor GAM42a_C 


Gammaproteobacteria 


GCCTTCCCACTTCGTTT 


23S 


1,027-1,043 




(26) 


Competitor BET42a_C 


Betaproteobacteria 


GCCTTCCCACATCGTTT 


23S 


1,027-1,043 




(26) 
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Table 3. Total cell counts and culturable population from the 
holothurian tissues 



Sample code a 


Total cell counts b 


Culturable cell 
counts 0 


PprppntatJP nf 

1 ^ [ \. >w 1 1 III \J l 

culturable cells (%) 


Si-1 


4.1xl0 6 


7.6xl0 4 


1.9 


Si-2 


6.1x10 s 


1.8xl0 4 


3.0 


Li-1 


3.6xl0 6 


8.8xl0 4 


2.4 


Li-2 


3.3xio 5 


2.6x10 s 


78.8 


U-l 


1.2xl0 6 


l.lxlO 4 


0.9 


U-2 


3.8x10 s 


1.7xl0 3 


0.5 


C-l 


9.9x10 s 


8.4xl0 3 


0.9 


C-2 


4.8x10 s 


1.3xl0 3 


0.3 



a Codes denote holothurian body parts and geographical origin. Si, 
Small intestine; Li, Large intestine; U, Surface; C, Coelomic fluid; 1, 
Funka Bay; 2, Ainuma Fishing Port. 

b Si-1, Si-2, Li-1, and Li-2, cells gf 1 ; U-l and U-2, cells cnr 2 ; C-l and 
C-2, cells mL 1 . 

c Si-1, Si-2, Li-1, and Li-2, CFU g" 1 ; U-l and U-2, CFU cnr 2 ; C-l and 
C-2, CFU mL 1 . 



100 



80 



60 



£ 40 

_ro 

o 



20 




Si-1 Si-2 Li-1 Li-2 C-1 C-2 U-1 U-2 

(n=29) (n=30) (n=27) (n=30) (n=29) (n=29) (n=28) (n=27) 

I I I I I I I I 

Small intestine Large intestine Coelomic fluid Body surface 



FJ Roseobacter clade 0 Shewanella spp. flU Pseudoalteromonas spp. 
Fj Vibrio spp. @ Flavobacteriaceae fj others 



Fig. 1. Composition of the culturable heterotrophs associated with 
holothurians. See Table 3 for sample codes. 



Culturable heterotrophs 

There were marked differences in the composition of 
culturable microbial community in each sample (Fig. 1). 
Heterotrophic microorganisms associated with holothurians 
were represented by members of eight different phylogenetic 



groups within the phyla Proteobacteria and Bacteroidetes: 
Roseobacter clade, Shewanella spp., Pseudoalteromonas 
spp., Vibrio spp., and family Flavobacteriaceae. The most 
frequently recovered population consisted of members of the 
genus Pseudoalteromonas (Fig. 1 and Table SI), which 
accounted for up to 52% of total culturable heterotrophs. The 
physiological tests of culturable populations indicated that 
bacteria with the ability to utilize alginate are more frequently 
recovered from digestive tracts than from other body parts 
(Table SI). Abundant and diverse heterotrophs were recov- 
ered from holothurian coelomic fluids. In the coelomic fluid 
of Funka Bay individual (C-l), Flavobacteriaceae members 
comprised 27.6% of isolates. In contrast, Vibrio members 
represented the second most frequently recovered population 
in the Ainuma coelomic fluid (C-2). Considering the total 
cell counts and the composition of culturable heterotrophs, 
microbial communities found in the coelomic fluids cannot 
be regarded as contaminants from other body parts. 

Microbial communities evaluated by 16S rRNA gene library 

Bacterial 16S rRNA gene clone libraries were successfully 
constructed using two universal primer sets from all holothu- 
rian samples (Table 1). The archaeal 16S rRNA gene was 
not amplified from any samples used in this study, although 
archaeal diversity was previously assessed for the midgut 
contents of a deep-sea holothurian species (27). A total of 
90 different bacterial phylotypes were identified from the 8 
libraries on the basis of classification with >94% identity 
(Table S2). The coverage values were 58.3 and 70.0% (small 
intestine), 61.5 and 65.2% (large intestine), 86.4 and 95.7% 
(coelomic fluid), and 61.5 and 86.4% (body surface). A 
distinctive bacterial community was detected in each holothu- 
rian niche (Fig. 2). 

The rRNA gene clones affiliated to the class 
Alphaproteobacteria were dominantly detected in all holo- 
thurian samples (13.0-54.5% in clonal frequencies). 
Alphaproteobacterial clones mainly belonged to three differ- 
ent subgroups (Fig. 3A). Roseobacter clade was mainly 
detected in small and large intestines (up to 43.5% in clonal 
frequency). These clones were closely related to clones or 
isolates previously retrieved from various marine environ- 
ments including coastal and pelagic seawater, sediments, and 
algae- and invertebrates-associated habitats (4, 5) (Table S2). 
Members of the genus Defluviicoccus were mainly detected 
from the body surface and the coelomic fluid (up to 54.5% 
in clonal frequency) (Table S2). Many alphaproteobacterial 
sequences from the coelomic fluid (28 and 37.5% in clonal 
frequencies; approximately half of the alphaproteobacterial 



Table 4. FISH-direct count analysis of metabolically active microbial populations associated with holothurians 



Samples 

Small intestine 
Large intestine 
Surface 
Coelomic fluid 



Total cell 
counts±SD 
(10 6 cells unir 1 ) 

22.1±3.0 
23.9±0.9 

2.2±0.3 

0.4+0.1 



EUB338 ARCH915 ALF968 



72.6±9.2 
74.419.8 
77.9±6.5 
76.314.4 



ND 
ND 
ND 
ND 



% of DAPI stained cellsiSD 

BET42a GAM42a EP402-423 CFB560 



27.813.9 
13.912.1 
46.917.5 
45.212.2 



7.010.8 
3.612.2 
3.913.5 
ND 



2.510.4 
3.411.4 
9.612.3 
2.510.3 



ND 
ND 
ND 
21.313.1 



4.011.2 
3.010.3 

ND 

ND 



PLA886 

12.210.2 
7.811.7 
ND 
ND 



SD, standard deviation («=4). 
ND, not detected. 
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Si-09 
(n=24) 



Si-40 
(11=23), 



Li-09 
(n=23) 



Li-40 
(n=26) 



U-09 
(n=22) 



U-40 
(11=23), 



C-40 

,(n=25) 



C-09 
(n=24) 



Small intestine 



Large intestine 



Body surface 



Coelomic fluid 



I Alphaproteobacteria ■ Deltaproteobacteria ■ Actinobacteria 
I Betaproteobacteria ■ Planctomycete | Firmicutes 
I Epsilonproteobacteria □ Fusobacteria □ Others 



Fig. 2. Similarity and composition of the bacterial population in 
holothurian tissues. The square distance (genetic similarity) was 
determined from the clonal frequency of each representative phylotype 
by the Ward method. Pie charts indicate the composition of bacterial 
population based on taxonomic grouping of 16S rRNA gene clone 
sequencing. See Table 3 for sample codes. 09: 24F-1509R primer set, 
40: 24F-1540R primer set. 



-Aquifex pyrophilus (M83548) 

-Defluvicoccus vanus (AF1 79678) 



I L 



strain Ellin5299 (AY234650) 
Oceanibaculum indicum (EU6561 1 3) 

clone B-09-325 (AB550448) 

-clone U-09-2 (AB550447) (+29 clones) 

E clone S-09-4 (AB560436) (+23 clones) 
-Ponticoccus litoralis (EF21 1829) 
r ropicibacter naphthalenivorans (AB302370) 

l clone B-09-239 (AB550452) 

T clone B-09-343 (AB550453) 

clone B-09-280 (AB550451) 
Rickettsia japonica (L36213) 
Orientia tsutsugamushi (U1 7256) 

clone B-09-288 (AB550460) (+12 clones) 



Rhodospirillales 



Rhodobacterales 



deep-sea clone BD2-8 (AB015538) 

— deep-sea clone EPR3967-02-Bc33 (EU491782) 
-Ehrlichia equi (M73223) 

Ehrlichia sennetsu (M7321 9) 

Wolbachia pipientis (X61 768) 



B 



-Desulfurella acetlvorans (X72768) 

- Wolineila succinogenes (M26636) 
-Helicobacter mustelae (M35048) 

-Helicobacter pylori (U01 330) 
— holothurian clone B-40-44 (AB550498) 
-cave clone PC061 1 0_B51_T7 (EU101204) 



rwolir 



I h, 



cave clone LKC3_22.72 (AY510167) 
Sulfuricurvum kujiehse (AB053951 ) 
r~ 'Candidatus Arcobacter sulfidicus' (AY035822) 

I {Arcobacter butzlerii (L1 4626) 

Arcobacter skirrowi (L1 4625) 



holothurian clone B-40-18 (AB550501) (+12 clones) 
deep-sea vent strain Am-N (AF357198) 
Sulfurospirillum arcachonense (Y1 1 561 ) 
-Sulfurospirlllum carboxydovorans (AY740528) 
Sulfurospirillum cavolei (AB246781 ) 
_x — Sulfurospirillum barnesii (AF038843) 
i- Sulfurospirillum deleyianum (Y1 3671 ) 

Fig. 3. Phylogenetic tree including representative holothurian clones 
as determined by neighbor-joining analysis. (A) and (B) were respec- 
tively constructed from 424 and 528 sites of the rRNA gene sequence 
that could be unambiguously aligned. Clones sequenced in this study 
are shown in red. The clonal frequency of each representative clone 
obtained in this study and DDBJ accession numbers are shown in 
parentheses. Branch points conserved with bootstrap value of >75% 
(solid circles) and with bootstrap values of 50 to 74% (open circles) are 
indicated. Some groups are represented by shaded trapezoids that 
indicate the numbers of sequences. Scale bars represent 0.05 sub- 
stitutions per nucleotide position. (A) Tree indicating the phylogenetic 
relationship among members of the Alphaproteobacteria. (B) Tree 
indicating the phylogenetic relationship among members of the 
Epsilonproteobacteria. 

clones) formed a novel cluster within the order Rickettsiales 
(Fig. 3A). Most Rickettsiales members have been recognized 
as obligate intracellular parasites of arthropods with the 
ability to infect vertebrates (44, 60). Recently, Rickettsiales 




Fig. 4. Epifluorescence micrograph of cells binding the Epsilonpro- 
teobacteria-specific probe (EP 402^123) in the holothurian coelomic 
fluid. Arrowheads indicate epsilonproteobacterial cells. Bar, 5 um. 



have also been found in the cells of leeches and marine 
ciliates (20, 60, 62). This study confirmed that Rickettsiales 
members could associate with more diverse animals than 
previously recognized. Clones affiliated to Planctomycetes 
and Fusobacteria were dominantly detected from holothurian 
digestive tracts (Fig. 2). These bacteria are frequently 
retrieved from marine snow and sediments (15), and thus 
may be transported with sediment particles into holothurian 
digestive tracts. Clones of the class Epsilonproteobacteria 
were detected only from the coelomic fluid (8.3 and 52.0% 
in clonal frequencies) (Fig. 2 and Table S2). As in the case 
of Rickettsiales clones, epsilonproteobacterial clones formed 
a novel clade which distantly related to members of the genus 
Sulfurospirillum (up to 92% sequence similarity) (Fig. 3B). 
Members of the genus Sulfurospirillum are microaerobic 
sulfur reducers, which include two validly described marine 
species, i.e., S. carboxydovorans and S. arcachonense (11, 
16). Known marine habitats from which Sulfurospirillum 
relatives were frequently recovered are restricted to sulfidic 
environments such as deep-sea vents and whale falls (32). 
Although the symbiotic ability of the Sulfurospirillum 
members is unclear, deep-sea vent Epsilonproteobacteria 
members have the ability to associate with various inverte- 
brates (30, 31, 33, 61). 

Whole-cell FISH analysis 

Microorganisms with high fluorescence emission after 
FISH hybridization have high growth rates and high metabolic 
activity (3). Metabolically active microbial populations were 
quantified using eight different fluorescent probes (Table 2). 
In general, FISH analysis confirmed the clone library data. 
Bacteria detected by the EUB338 probe accounted for 
approximately three fourths of total cells in all holothurian 
samples (Table 4). Epsilonproteobacteria specific probe- 
binding cells were observed only in the coelomic fluid (Fig. 
4). ALF968 probe-binding cells were abundantly found in 
all holothurian samples. These cells existed freely even in 
the coelomic fluid where Rickettsiales members accounted 
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for half of the alphaproteobacterial clones, suggesting that 
Rickettsiales members detected in the coelomic fluid lack the 
ability to inhabit holothurian cells. 

Conclusion 

This study revealed for the first time the diverse and 
abundant microbial community associated with the holo- 
thurian body. In particular, it was remarkable that probably 
metabolically active and phylogenetically unique Epsilon- 
proteobacteria and Rickettsiales members were discovered 
in the holothurian coelomic fluid, which contains antibacterial 
compounds (8). In addition, Echinodermata immunocytes 
produce reactive oxygen species (ROS) when in contact with 
bacteria or bacterial cell wall proteins (28). Nevertheless, 
potentially pathogenic microbes and fungi were repeatedly 
isolated from the holothurian coelomic fluid (1, 46). In 
addition, the coelomic fluids of oysters and abalones were 
used as inocula for isolating novel microorganisms (10, 24, 
52). Echinoderms, at least those inhabiting intertidal areas 
where physical and chemical conditions fluctuate, have 
increased ability to maintain homeostasis in coelomic fluid 
(6). These suggest that coelomic fluid of marine invertebrates 
provides various oceanic microorganisms with unique, 
stable habitats. Focal points raised by this study for future 
research include: time course and evolution of symbiont 
diversity, microbial metabolism, host-microbe interactions 
and potential biotechnological and fisheries implications. 
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